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Cosmic Rays 


High-energy protons and heavy 
ions 

Very broad range of energies 

Galactic and Extragalactic 

maximum energy predicted at 
60x1018 eV (= 60 EeV) (GZK- 
cutoff) 
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Extreme Energy Cosmic Rays 

(EECRs) 
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• EECRs (E > GOEeV) suffer almost no deflection from 
magnetic field in extraga lactic, galactic, and solar system 

• They point back directly to the location of their original 
sources "Charged Particle Astronomy" 


• Detection of EECRs gives us information about their 


sources 
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Question 6: 

How Do Cosmic Accelerators 
Work and What Are They 

Accelerating? 
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^''Connecting Quarks with the Cosmos: Eieven Science Questions 
for the New Century Committee on the Physics of the Universe", 
National Research Council , ISBN: 0-309-50569-0 (2003) 
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Spectra of leading EECRs observations 

• EECRs Flux is low of the order of 1 particle/kmVsr/century 

* At high end of the spectrum reduces to 1 particle/kmVsr/millennium! 
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Extreme Universe Space Observatory 
Mission 




• It will be the first space- 
based observatory to 
use the Earth's 
atmosphere to discover 
the origin of EECRs. 
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• A collaboration of 15 
countries, 85 
institutions and ~ 333 
scientists 
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• Studies of the nature of EECRs are based on 
the measurement of EASs, 

• EASs are cascades of secondary particles in 
the atmosphere as a result of the interaction 
of EECRs with the Earth's atmosphere 


* Maximum atmospheric 
depth depends on: 
primary energy, 
nature of the primary 
particle 

details of the interactions 
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• Atmosphere is modeled as cylinder of radius 0.5 km 
and divided into 112 layers of 1 km thickness and 
varying density 
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EEL Modeling EAS using Geant4 


• Using Geant4 model FTF_BIC 

Proton and Iron as primaries 

• Time and data size are huge! 

• Simulating 1 event of ITeV proton: 

> Time ~ 20 hrs 

> Size -200GB 

• Simulation time for EAS scales 
with energy of primary particle 
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• Let be energy of primary particle, and E be 
energy of secondary particle 

Thinning level 

• All secondaries with energies greater than EtH 
are followed, BUT if the energy sum of all J 
secondaries produced in a certain Interaction 
falls below thinning energy 


^th^o 



only one particle is followed 
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Geant4 vs CORSIKA 
lOOTeV Proton 
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Geant4 vs CORSIKA 
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• Thinning algorithm cuts simulation time 
significantly 

• X^Q^is not affected by thinning level 

• Location of depends on primary particle 


• Preliminary calculations show that Geant4 
predicts low X^^^ compared to CORSIKA: 

> Magnetic Field 

> Geometry 


4/7/2015 


18 



NASA 


POSTDOCTORAL PROGRAM 

Administered by Oak Ridge Associated Universities 


Acknowledgement 



Marshall 

Space 

Flight 

Center 


This work is supported by an appointment to the 
NASA Postdoctoral Program at Marshall Space 
Flight Center, administered by Oak Ridge 
Associated Universities through a contract with 


NASA 



www.orau.org 


4/7/2015 


19 





